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ABSTRACT 

We present results from a spectre- interferometric study of the Miras o Cet, R Leo and W Hya 
obtained with the Keck Aperture Masking Experiment from 1998 Sep to 2002 Jul. The spectrally 
dispersed visibility data permit fitting with circularly symmetric brightness profiles such as a simple 
uniform disk. The stellar angular diameter obtained over up to ~ 450 spectral channels spaning the 
region 1.1-3.8/im is presented. Use of a simple uniform disk brightness model facilitates comparison 
between epochs and with existing data and theoretical models. Strong size variations with wavelength 
were recorded for all stars, probing zones of II2O, CO, OH, and dust formation. Comparison with 
contemporaneous spectra extracted from our data show a strong anti-correlation between the observed 
angular diameter and flux. These variations consolidate the notion of a complex stellar atmosphere 
consisting of molecular shells with time-dependent densities and temperatures. Our findings are 
compared with existing data and pulsation models. The models were found to reproduce the functional 
form of the wavelength vs. angular diameter curve well, although some departures are noted in the 
2.8-3.5 /im range. 

Subject headings: instrumentation: interferometers - techniques: interferometric - stars: AGB and 
post-AGB - stars: individual: Mira stars: individual: W Hya stars: individual: R 
Leo 



1. INTRODUCTION 

Due to their large diameters and high luminosities, 
Mira variables are a favourite target for optical and near- 
infrared (NIR) interferometric observations, which have 
shown dependencies of diameter on wavelengt h, pul- 
satio n phase and puls at ion c ycl e (e.g., iHaniff et al.l 



19951 ' Ivan Belle et all 119961:'' iPerrin et al.l Il999t 
Jacob et al.l l2000t lYoung et all 120001: iHofmann et al l 



20021: iThomps on et a ll 120021: [Ireland et al.l l2004d : 
Woodruff et alT 2004: W oodruff et al.ll2008D . This wealth 



of interferometric information has advanced the studies 
of m olecular and dust a bu ndances in the atmosphere 
(e.g. Ilreland et al]l2004al lbl: llreland fc Scholz 200 61), the 
pulsation mo d e of t hese stars (e.g.. I Woodruff et al.ll2004l : 
iFedele et al.l 120051), and photqspher ic / circumstellar 
asymmetries (e.g. iRagland et al.l l2006( ) as well as the 
characteristics of the circumstellar environment (e.g. 
iDanchi et a"l]|1994D . 

Miras are long period variables (LPVs) of roughly 
one solar mass which evolve along the Asymptotic Gi- 
ant Branch (AGB) and are characterized by well-defined 
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pulsation periods and large pulsation amplitudes (up to 
AV^ K, 9). The extended atmospheres of these late- type 
giants show phase-dependent density and temperature 
stratifications which, together with complicated inten- 
sity profiles and the presence of molecular layers, leave 
no simple observable quantity to define the structure 
of the star's envelope. The pulsations also help drive 
strong winds thus enriching the interstellar m edium with 
mole c ules such as H2O, CO, T iO and SiO (Tsuii et all 
119971: iTei et al.ll2003at lOhnakal [2004) . A better under- 
standing of their atmospheric structures will thus further 
the understanding of these late stages of stellar evolution, 
and might help shed light on one of the most important 
sources for the chemical enrichment of the interstellar 
medium. 

To fully characterize the structure of a Mira star's 
atmosphere, one would need the complete intensity 
profile at all wavelengths and all pulsation phases, 
recorded over a large number of pulsation cycles. The 
last decade has shown a multitude of angular diameter 
measurements of Miras in various wavelengths, with 
some instruments returning imp ressive resolution of 
the brightness prof ile (e. g IOT A: [Ragland et al.ll2006l . 
VLTI: Wittkowski et al.1 l2008l etc). The difficulty 
lies in gathering multi-phase, multi-wavelength data 
on a single object. Here, we present multi-epoch, 
spectro-interferometric observations of the nearby 
Mira stars o Cet, R Leo and W Hya covering the 
NIR spectrum from 1.1-3.8 /xm near-simultaneously, 
with one-dimensional spatial information out to the 
diffraction limit of the Keck I telescope. With our data 
set we are able to sample phase-dependent molecular 
stratification and help probe the structure of the stellar 
atmospheres. 
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TABLE 1 
Observed objects 



Name 


Period 


Spectral 


Distance 




[days] 


Type Range 


[pc] 


W Hya 


385 


M7.5-9ep 


104 ± 12 


R Leo 


312 


M6-9.5e 


111 ± 17 


o Get 


332 


M5-9e 


92± 10 



References. — Object Period from the 
American Association of Variable Star Ob- 
servers (AAVSO) visual light curves (A. A. 
Henden et al. 2006, private communication), 
M Spe ctral Type Range from Sloan & Pric3 
1 1998 ) and distances from Ivan LeeuwenI 
1 20071 ) 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. Aperture-masking Observations 

Our sample of 3 Miras (W Hya is sometimes classified 
as a semi-regular pulsator with strong Mira characteris- 
tics) were chosen for their large angular diameters and 
NIR brightness (see Table [Ij . Observations were per- 
formed with the 10 m Keck I telescope using the Near 
Infrared Camera (NIRC). The telescope pupil was con- 
verted into a sparse interferometric array by placing an 
aperture mask in the beam in front of the infrared sec- 
ondary mirror, allowing the recovery of the Fourier am- 
plitudes and closure phases for baselines up to 10 m. 

Our experimental methods are derived from the highly 
successful aperture masking program at the Keck obser- 
vatory which, in ad dition to studies of evolved pulsat- 
ing s t ars and giants (Tuthill et al.lf2000al : iMonnier et al.l 



[200l iWoodruff et al 



2008i, has also delivered advances 



in stellar astrophysics ranging from youn g stellar objects 
(|Danchi et al.ll200lHTuthill et al.ll2002bl ) to dusty mass- 
loss shrouds in proto-planetary nebulae and mass ive stars 
(|Monnier et al.l[2000l : Irilthill et al.ll200"2aL [200l . A full 
description of the experiment, including a discussion of 
the conceptual principles and signal-to-noise consider- 
ations under l ying m asking interferometry, is given in 
iTuthill et al.l (|2000bD . while further discussion of sys- 
tematics and seeing induced errors can be found in 
IMonnier et all (|2004f ). 

In a significant extension to the capabilities of our pre- 
vious experimental setup, data presented here were ob- 
tained in a spectrally-dispersed manner delivering con- 
tinuous wavelength coverage across the infrared J, H 
and K bands (from 1.0 /im to 3.7 /im). The use of a 
one-dimensional non-rcduntant aperture mask, as de- 
picted in Figure[Tl in combination with a cross-dispersing 
;rism element already installed in the NIRC camera 
Matthews et al.|[l996r ). resulted in spectral information 
being encoded in one dimension on the readout array, 
while spatial interference fringes were recorded in the 
other dimension. Note that with this setup, the full two- 
dimensional imaging data recovered in the previous ex- 
periment was sacrificed for one-dimensional brightness 
profiles across the waveband. 

Figure [1] shows the aperture mask, which was designed 
to fit within the Keek's segmented primary. The design is 
a trade-off between uniform Fourier coverage and avoid- 
ing features within the pupil, such as the central obstruc- 



tion and the boundaries between mirror segments. Al- 
though the mask only passes a small amount of the light 
gathered by the 10 m primary mirror, and still more light 
is lost at the slit entrance, the LPV targets discussed in 
this paper are among the most luminous objects in the 
near-infrared sky and the low overall throughput did not 
adversely affect the signal-to-noise of the final data. On 
the contrary, we found that there was even sufficient flux 
recorded on the shoulders and in between the traditional 
near- infrared band windows, where atmospheric absorp- 
tion is high, so that nearly continuous wavelength cover- 
age could be obtained, with the exception of the strong 
telluric emission features, due to water, between 2.6 and 
2.9 Aim. 

The NIRC camera has three grisms intended to cover 
the J-H, H-K and K-L spectral regions (see Table [2]). All 
three of these grisms were used and, as the reasonable 
degree of wavelength overlap between them provided a 
useful check on the reliability of the data. Also, four dif- 
ferent slit widths (1.5, 2.5, 3.5, 4.5 pixels as projected on 
the detector) were used in order to strike a balance be- 
tween spectral resolution and available flux. The grisms' 
wavelength resolution are given in Table [2] and are to be 
multiplied with the slit width used in the observation (see 
Table [3]) convolved with 1 pixel of the detector to obtain 
the "true" wavelength resolution of the instrument. 

After passing the mask and being dispersed in the 
NIRC spcctrographic mode, the beams were focused on 
the 256 X 256 pixel array at a scale of 20.57 mas pixel"^. 
This is sufficient to Nyquist-sample data collected for 
2 /im or longer wavelengths. An example of a typical 
short-exposure (ri„t=0.14s) frame is sho wn in the left 
panel o f Figure El Following the methods of l Tut hill et al.l 
()2000bf l. 100 of these frames could be collected with rea- 
sonable efficiency (duty cycle ^^20%) into a data cube for 
processing, to statistically calibrate the atmospherically 
degraded point-spread function. 

The power spectrum accumulated over a cube of 100 
frames is illustrated in the right panel of Figure [2l Power 
on fifteen baselines can be clearly seen as peaks appearing 
at discrete spatial frequencies that vary smoothly with 
wavelength. By recording the power on each baseline 
for a given wavelength bin, a one-dimensional visibility 
function could be recovered for that channel. When cal- 
ibrated for the total flux level in each channel, these raw 
1-dimenional V"^ data could be collated and fitted with 
any convenient model (e.g. a uniform disk; model fit- 
ting is discussed in Section r2.2|l . yielding a continuous 
sampling of the angular size of the target as a function 
of wavelength. Before such fitting could take place, cal- 
ibration for the telescope-atmosphere transfer function 
was needed. For a number of instrumental reasons, this 
procedure turned out to be more complicated than in the 
past, as described below. 

At wavelengths shorter than 2 /im, the longest base- 
line corresponds to spatial frequencies higher than the 
Nyquist sampling condition. This power is aliased back 
into the power spectrum at lower spatial frequencies. 
This effect can be observed in the power spectrum in Fig. 
[2] as an apparent reflection of the long-baseline power 
peaks when the wavelength is below w 2 /im. Where this 
wrapped signal overlaps with power from shorter base- 
lines, the data become confused and were discarded. 

In common with established practice in interferometry, 
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TABLE 2 

Table of gmsms 



Grism Keck/NIRC blocking 


Grism wavelengths 


effective wavelengths 


spectral resolution 


Filter Name 


[^tm] 




[10-'-^ fim] 



GR150 JH 0.82-2.08 1.0-1.6 ~5 

GR120 HK 1.05-2.58 1.4-2.5 ~6 

GR60 KL 1.98-4.75 2.1-4.5 ~ 10 



TABLE 3 

Observations of Mira stars in our sample, including cycle and visual 
phase $ (cycle=0 for first observation, visual phase=0 at maximum 

light) 



Object name 


Date 


JD 


-I- 


Grisms 


Slit width 




-2450000 






[pixels] 



o Get 


1998Sep29 


1056 


0.71 


GR60,120,150 


1.5, 1.5, 1.5 




2002Jul23 


2479 


4.98 


GR60,120 


1.5, 3.5 


R Leo 


1999Feb04 


1213 


0.49 


GR120,150 


3.5, 3.5 




1999Apr25 


1295 


0.75 


GR60,120,150 


4.5, 4.5, 4.5 


W Hya 


1999Feb05 


1213 


0.58 


GR60,120,150 


2.5, 3.5, 3.5 




1999Apr25 


1295 


0.79 


GR60,120,150 


4.5, 4.5, 4.5 




2000Jan25 


1570 


1.53 


GR60,120,150 


1.5, 3.5, 4.5 



TABLE 4 

Reference stars with spectral types and estimated 
diameters 



Reference star 


Spectral 


Adopted UD 


Reference 




Type 


Angular Diameter [mas] 




2 Gen 


M4.5HI 


14.7 


1 


a Get 


M1.5HI 


11.7±0.6 


2 


a Hya 


K3ii-ni 


9.1 ±0.1 


3 


a Lyn 


K7III 


7.2 ±0.6 


2 


a Tau 


K5HI 


19.7 ±0.1 


4 



Ref ere nces. — (1) |Du mm & Schild ( 199 ^); (2 
119981) : (3) IMozurkewicheTaLi m03, ): (4) ,Perrin et al.l (119981 ) 



get and hence permits fitting with circularly symmetric 
brighness profiles such as a simple uniform disk. The 
spectro-interferometric data were recorded at 6 different 
epochs spanning more than 4 years. Table [3] lists these 
observations. 




observations of PSF calibrator stars with smaller and 
well-characterized angular diameters were interleaved 
with those of the science target (see Table 2]). However, 
for observations with the grism setup utilized here, an 
observational difficulty compromised this "standard" cal- 
ibration process. Due to the experimental requirement 
for a stationary optical path between the mask and detec- 
tor, the Keck image rotator was not used for any masking 
experiment. This precluded use of the normal telescope 
guiding system, and each series of 100 exposures (last- 
ing of up to a couple of minutes) were taken with the 
telescope in a blind tracking mode. While small drifts in 
tracki ng had no impact o n the original masking experi- 
ment (jTuthill et al.|[2000b( ). this was not the case when 
the grism was employed because wander of the stellar 
image across the slit was found to result in significant 
modification to the optical transfer function. This moti- 
vated the construction of an alternate calibration strat- 
egy outlined in section 12.21 

The spectrally dispersed visibility data resolve only one 
spatial dimension at a single position angle for each tar- 



FlG. 1. — Pupil geometry used for the Keck masking + grism 
experiment. A 6-hole non-redundant linear array (dark squares) is 
shown overlaying a scaled version of the segmented Keck primary 
mirror. The boundary of the f/25 secondary mirror as projected 
on the primary is represented as the bold black line. 



2.2. Uniform Disk Diameters and Calibration 

As discussed in [Woodruff et al.l (|2008D . the true stellar 
intensity profile is not a uniform disk (UD), but shows a 
complex center-to-limb variation that varies with wave- 
length, pulsation phase and cycle. Fitting the observed 
1-d visibilities with a simple profile like a UD provides 
a useful estimate of the apparent change of size as a 
function of wavelength. Since we are only resolving low- 
resolution structure in the star's intensity profile, it be- 
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pixels Spatial Frequency <rad") 

Fig. 2. — Spectrally dispersed interferogram of W Hya taken with the GR120 grism (left). The interference pattern generated by the 
6-hole linear mask (see Fig. [TJ along the abscissa is dispersed with wavelength along the ordinate. The power spectrum (right) is the 
result of averaging the power spectra of 100 such data frames. Note the regions of attenuated signal due telluric absorption in both the 
interferogram and the power spectrum. 



comes difficult to differentiate between a UD, a fully- 
darkened disk, or a Gaussian. Any one of these simple 
models will show a similar relative change in diameter 
with wavelength, so we chose UD diameters to allow com- 
parison of findings with existing literature and to avoid 
the difficulties en countered when fitting more soph isti- 
cated profiles (cf. iHofmann et ZI[T998t ISchoyl2003h . 

All spectrally-dispersed observations using the grism 
were accompanied on the same night with orthodox two- 
dimensional aperture masking observations in the NfRC 
narrowband filter set (with the exception of the L band 
for W Hya in 1999 Apr 25). Since calibration of these 
observations was not impacted by the blind tracking er- 
rors, we could obtain good contemporaneous angular di- 
ameters of our targets at a range of discrete wavelengths 
across the grism passband. Assuming that the tracking- 
induced instabilities in the optical transfer function were 
a smoothly-varying function with wavelength, we could 
use our robust discrete UD angular diameters to pin the 
calibration at those wavelengths, and perform a smooth 
interpolation in between to cover all passbands used in 
the dispersed experiment. In the case of W Hya in 1999 
Apr 25, we used L3.08 2-D making observations con- 
ducted at t he same phase albeit o ne year earlier, on 1998 
Apr 14 (cf. iWoodruff et"an[200l . for one calibration of 
the GR60 grism. 

Observations of the PSF reference stars were then used 
as check stars to verify the utility of this calibration pro- 
cess. Although these check stars were typically smaller 
and fainter in the NIR than the targets, and therefore 
yield proportionately higher relative errors on the di- 
ameter measurement, the results of applying our cal- 
ibration methodology to these objects confirmed that 
it was robust. PSF reference stars had typical spec- 
tral types from late K to early M, and (in comparison 
with Mira stars) these stars have compact atmospheres 
so that any diameter changes with wavelength are ex- 
pected to be small. After applying our new calibration 
strategy to these check stars, the results verified these 
expectations and showed that the method produced a 
well-characterized constant diameter measurement over 
all three grism passbands. This is shown in Fig. [3]for the 
reference star 2 Gen, observed contemporaneously with 



W Hya in 1999 Apr 25. 

2.3. Low resolution spectra 

The observations of reference stars also allowed the 
spectral response of the optical system to be estimated. 
We extracted the total flux as a function of wavelength 
for both science objects and reference stars, making it 
possible to recover basic low resolution NIR spectra in 
the 1.0-4.0 /im range (see Table[2]) contemporaneously to 
our visibility measurements. Since the spectral calibra- 
tion suffers from the same observational difficulty which 
compromised the calibration process of the visibilities 
outlined in Section [21 it was not possible to obtain cali- 
brated absolute fluxes from our measurements. 

3. RESULTS 

Figures m to [9] show the UD diameters of W Hya, R Leo 
and o Get as a function of wavelength (UD(A)), together 
with contemporaneous spectra. The gaps in the data 
correspond to regions of attenuated signal due to the tel- 
luric absorption bands, with the exception of the gaps 
at wavelengths greater than ^3.5/zm, which are caused 
by saturation of the detector due to high flux from the 
sky background. Also shown in the Figures are the best 
fitting UD(A) from dynamic model atmospheres with the 
respective model spectra which are described and com- 
pared with the data in Section IH 

The UDs show wavelength-dependent features that are 
consistent with expected opacity changes due to molec- 
ular abundances. Although some VO is responsible for 
UD(A) variations in the J-band, the majority of the fea- 
tures are caused by the presence of H2O (plus contribu- 
tions of CO and OH) in the stellar atmo sphere (see, e.g., 
iTei et al.l[2003bHWittkowski et al.ll2008D . There is a very 
close anti-correlation between the spectral features and 
the apparent angular diameters at varying wavelengths. 
We find local minima in the UD diameters at ~ 1.3-1.4, 
1.6-1.7, and 2.2 /im, corresponding to the regions of in- 
creased flux in the spectra. This anti-correlation can 
be explained by a molecular blanketing model, where 
opaque layers above the continuum-forming layers (the 
photosphere) lead to a perceived UD angular diame- 
ter increase at certain wavelengths. This wavelength- 
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dependent absorption is also responsible for the dimin- 
ished flux in the same bandpasses. The less contaminated 
bandpasses, which are closer to sampling the continuum- 
forming layers, exhibit smaller measured UD diameters 
and higher fluxes. These results can be readil y com- 
pared with the findings for S Ori of iWittkowski et al.l 
|2008), which show UD diameter minima in the same 
wavelength regions and present a very close match re- 
garding the overall form of the UD(A) curve. 

For the partially-resolved reference s tar 2 Cen we as- 
sume d a UD diameter of 14.7 mas (cf. iDumm &: Schildl 
Il998f l at the same wavelengths used to pin the Mira grism 
data (cf. Figure [3|). Note that some variation of size with 
wavelength is expected for this M4.5 III star. The corre- 
sponding spectrum is also fairly smooth, exhibiting only 
small H2O absorption features at ^1.4 and 1.9 /im. 

Figures [Hand [7] also show significant cycle-to-cycle and 
phase effects, which tend to be more pronounced towards 
shorter wavelengths. This effect is predicted by the mod- 
els, as discussed in the next Section. 

4. COMPARISON WITH MODELS 

In this section the predictions of recent dynamic model 
atmospheres of M-type Miras are compared with the 
present observations. The distances to the model stars 
were chosen such that the UD(A)inodoi overlaps our grism 
data in the wavelength region 2.2-2.5 /j,m. These dis- 
tances do not guarantee good fitting to model photom- 
etry since we only consider a comparison between the 
functional form of the model predictions as a function 
of wavelength and the observed data. A comparison be- 
tween model photometric and angular diameter fit dis- 
tances will be given in a future paper describing the mod- 
els. The new models, spanning several pulsation periods, 
are based on a new code for treating pulsation and cal- 
culating atmos pheric stratifica t ions, as described below. 

The code of iKeller fc modi (pOfffih was used for con- 
structing two series of self-excited pulsation models with 
parameters close to the well-studied prototypes o Get and 
R Leo: (i) o54 model series with period ~ 330 days (non- 
pulsating "parent star" with mass 1.1 Mq, luminosity 
5400 Lq, radius 216 Rq, effective temperature 3370 K, 
solar abundances with Z=0.02); (ii) r52 model series with 
period - 307 days (1.1 Mq, 5200 Lq, 209 Rq, 3400 K, so- 
lar Z=0.02). The temperature structure of the atmo- 
spheric layers was calculated with an opacity-sampling 
technique accounting for major molecular species blan- 
keting M-type atm ospheres as well as for dust opacity 
(|Ireland et al.l[2008L henceforth called ISW08). The o54 
series covers three time intervals, each containing 1.25 
(interval hx, hereafter o54(hx)), 1.25 (jx) and 3.75 (fx) 
successive cycles in phase steps of about 0.1 times the 
pulsation period. For the r52 series, there are two time 
intervals with 2.25 (gx) and 2.5 (fx) successive cycles. 

Instantaneous relaxation of hot matter behind shock 
fronts is assumed, and the approximation of local ther- 
modynamical equilibrium was adopted except for scatter- 
ing processes affecting continuous, TiO line and dust ex- 
tinction (see ISW08 for details). All model atmospheres 
are cut off at 5 parent star radii, which is considered to 
be the transition zone to a wind-dominated circumstellar 
shell. Though W Hya has a longer period than these pro- 
totype model series, we would expect that scaling their 
atmospheric geometry structure to the larger size of a 



longer-period Mira provides a fair approximation for this 
star. 

Inspection of the model series shows that, apart from 
the strong dependence of the shape of UD(A)i„odoi curves 
with optical phase, these curves also vary quite notice- 
ably between cycles. Whilst an observed small set of iso- 
lated UD(A)obs values or an observed narrow wavelength 
interval of the UD(A)obs curve might readily be fitted by 
model predictions in different cycles, close agreement is 
much harder to achieve over a broad wavelength range, 
as we find with our NIR observations. This reflects the 
varying structure of upper atmospheric layers, depending 
on the shock front position in the present cycle as well 
as on the shock front history of preceding cycles, which 
results in varying strengths of different molecular absorp- 
tion features. Systematic observations of monochromatic 
cycle-to-cycle size variations of Miras are not available, 
but the present observations of W Hya in 1999Feb (phase 
0.58) and 2000Jan (phase 1.53) are almost one full cy- 
cle apart and clearly show noticeable differences of the 
UD(A)obs curve, which are significantly larger than the 
small 0.05 phase-difference effect expected from models 

(Fig. a. 

4.1. W Hya 

We compared the 3 phases of observation of W Hya 
with predictions of different cycles of both the o54 and 
r52 model series. Whilst no good agreement was found 
with available r52 cycles, there is reasonable agreement 
with one modelled phase sequence of the o54 series when 
the phase assignment of the model series is estimated to 
be slightly earlier, say of the order of 0.05 to 0.1 (which is 
well within the typical model vs. star phase uncertainty 
of at least 0.1). Even better agreement is achieved be- 
tween our phase 0.58 observation and models with nearby 
phases 4.55 and 4.61 but, since subsequent model cycles 
are presently not available, the observed phases 0.79 and 
1.53 cannot be compared to models in this case. 

For good agreement with observed UD diameters, the 
model star had to be put at a "fit-distance" of about 
78 pc, as compared t o 104 ± 10 pc for the real star 
(jWhitelock et al.ll2008f ). where the difference reflects the 
fact that W Hya has a noticeably longer period (385 days) 
and hence is larger than the o54 model star (330 days, 
parent star radius 216Rq). 

In Figure H we show the 1999Feb and 2000Jan ($ 
0.58 and <!> = 1.53, respectively) observation together 
with UD(A) curves of o54(fx) models with phases 2.51 
and 3.49 of successive cycles. Figure[5]shows the 1999Apr 
observation ($ — 0.79) together with the phase 2.70 
model, as well as the phase 2.80 and 3.70 models for 
demonstrating significant phase and cycle-to-cycle ef- 
fects. Figure [6] shows the 1999Feb observation of W Hya 
($ = 0.58) with phase 4.55 and 4.61 models plus the 
phase 3.59 model of the preceding cycle which gives a less 
satisfactory overall fit. We note that the cycle contain- 
ing phases 3.49 to 3.70 has generally smaller UD values 
than the other cycles used in Figs. [3]to[Sl i.e. the model 
star appears smaller in this cycle and had to be scaled 
accordingly. 

4.2. R Leo 

For the observations of R Leo at two phases, none of 
the four available cycles of the r52 model series (which 
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GR60 O 

2 Cen 1999Apr GR150 • 

2.2 /xm UD diameter 




1.5 



2.0 2.5 3.0 

Wavelength [/J^m] 



3.5 



4.0 



Fig. 3. — Top: UD angular diameters for the reference star 2 Cen as a function of wavelength observed on 1999 Apr 25 with the 
GR150, GR1 20 and GR60 grisms f blue, green and red diamonds respectevely ) . Th e magenta filled squares show the 2.2 /^m diameter of 
14.7 mas (cf. IDumm Hi Schildlll99 g) which are used for calibration (see Section 12.21 ) Bottom: Low resolution spectra for 2 Cen observed 
contemporaneously with the GR150, GR120 and GR60 grisms (blue, green and red lines respectively). The flux per unit wavelength is in 
arbitrary units and normalized at the H-hand. 



was originally designed to describe R Leo) matches the 
observed UD(A) curve well, whereas two cycles of the 
o54 series with a slightly longer period give a fair fit. 
For good agreement with observations, the model star 
was set at a distan ce of ^ 93 pc, to be compared with 
parallaxes given by iGatewoodl ()1992l 120 ± 14 pc) and 
Ivan LeeuweiJ ([2007L HIPPARCOS 73 ± 6pc) and with 
a weighted-mean distance suggested by iWhitelock et alj 
(|lb08. 111 ± 17pc). Here, similar to the case of W Hya, 
we note that the $ — 0.75 observation is in the phase 
range of pronounced UD(A) phase dependence when the 
Mira light curve increases very steeply. Fig. [7] shows 
the observed UD(A) curves together with model curves. 
These give a fair fit within the HKL bands, whereas the 
modeled UD diameters in the J bandpass is about 15 
percent larger than the observed value at phase 0.49. 

4.3. o Get 

For o Get, observations were made at two phases in 
two quite distinct cycles (1998Sep and 2002Jul) so that 
the phase-cycle-effects cannot be studied readily in this 
case. Two basic problems occur when we fitted available 
models to the observed UD(A) curves, (i) The model par- 
ent stars have to be placed at a distance of ~80 pc (o54) 
or ~70pc (r52), both of which are marginally smaller 
th an values given in the literature, e.g. 92 ±10 pc af- 
ter [WhiteiockiraD (|2008'). (ii) The decrease of the UD 
diameter measured on 1998Sep at the long-wavelength 
side of the strong 3 /im water band tends to be notice- 



ably steeper than predicted by any models. Figure [8] 
shows the 1998Sep observation of o Get at $ = 0.71 to- 
gether with same phase models of two successive cycles, 
showing cycle effects, and a model with 0.1 later phase, 
demonstrating the strong phase effects around $ = 0.8. 
In Fig. [5] the (<i> = 4.98) observation of o Get is shown 
with two model curves of the r52(fx) series. 

5. SUMMARY 

We have measured the wavelength dependent diam- 
eters of three Mira variables stars from 1.0-3. 4 /im and 
their respective spectra at up to 3 different epochs, 
the first study of this type in the near-infrared. All 
observed stars show strong variations of their UD 
angular diameters as a function of wavelength, often 
exhibiting a factor of ~ 2 in UD diameter between 
1.0 /xm and 3.0 /im. We find that the UD(A) relation 
shows variations with phase and the pulsation cycle, 
revealing the complexity of cycle and phase-dependent 
opacity contributions from molecules, predominantly 
H2O and GO, in different layers. 

We found that there is fair agreement between the mea- 
sured UD(A) relationship and spectra, and those pre- 
dicted by theoretical models designed to represent o Get 
(see ISW08). When comparing the two stars whose pa- 
rameters are thought to resemble the model input param- 
eters the most. R Le o and o Get, with model predictions 
(jlreland et al.ll2008f ). we find that the model UD diam- 
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Fig. 4. — Top: UD angular diameters for W Hya as a function of wavelength observed on 1999 Feb 5 ($ = 0.58) and 2000 Jan 25 
(# = 1.53) with the GR150, GR120 and GR60 grisms (blue, green and red diamonds respectively for 1999Feb, and grey diamonds for 
Jan2000). The black solid and dashed lines show the fairly well fitting models of the o54(fx) series at phases 2.51 and 3.49. The black 
and magenta filled squ ares show the narrow-band two-dimensional UD angular diameters observed on the same nights which are used for 
calibration (see Section 12.21 ) Bottom: Low resolution spectra for W Hya observed contemporaneously with the GR150, GR120 and GR60 
grisms (blue, green and red lines respectively for 1999Feb, and grey lines for Jan2000). The flux per unit wavelength is in arbitrary units 
and normalized at the //-band. The black line shows the spectrum predicted by the model of the o54(fx) series at phase 2.51 (for the sake 
of clarity and since both model spectra differ only slightly, only one model spectrum is shown). 




Fig. 5. — Same as Fig. l4lbut for W Hya observed on 1999 Apr 25 ("3? = 0.79) only. Also shown is the UD(A) curve as predicted by the 
stellar models discussed in Section |4] as a black solid, dotted and dashed lines for the phases 2.70, 2.80 and 3.70 from the o54(fx) model 
series, respectively. The model spectrum corresponds to the phases 2.70 from the o54(fx) model series. 
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Fig. 6. — Same as Fig. |4]but for W Hya observed on 1999 Feb 5 ('I' = 0.58) only. Also shown are the UD(A) curves as predicted by 
the stellar models discussed in Section |4] as black solid, dotted and dashed lines for the phases 4.55, 4.61 from the o54(fx) model scries 
as well as the model at 3.59 of the preceding cycle for demonstrating significant cycle-to-cycle effects, respectively. The model spectrum 
corresponds to the phases 4.55 from the o54(fx) model series. 



60 



50 



a 40 



30 



0.5 





1999Feb 


GR120 


o 


1999Feb 


GK 1 50 




1999Apr 




■ 


1999Feb 


2-D 


■ 


1999Apr 
- models 


2-1) 




R Leo 

1999Keb and 1999Apr 




1.0 



1.5 



2.0 2.5 

Wavelength [/.tm] 



3.0 



3.5 



4.0 



Fig. 7.— Same as Fig. Hbut for R Leo observed on 1999 Feb 04 ($ = 0.49) and on 1999 Apr 25 (* = 0.75) with the GR150 and GR120 
grisms (blue and green symbols respectively for 1999Feb, and grey symbols for 1999Apr). The model UD(A) curves shown as a solid and a 
dashed black line correspond to the phases 2.51 and 2.70 of the o54(fx) model series. The model spectrum corresponds to the phases 2.70 
from the o54(fx) model series. 



Spectro-interferometry of Mira Variables from 1.1 to 3.8 microns 



9 




Fig. 8. — Same as Fig. |4]but for o Get observed on 1998 Sep 29 ($ = 0.71). The model UD(A) curves shown as soUd, dotted, and dashed 
black lines correspond to the phases 0.70, 1.70 and 1.80 from the o54(hx) model series, respectively. The model spectrum corresponds to 
the phases 0.70 from the o54(fx) model series. 




Fig. 9. — Same as Fig. |4]but for o Cot observed on 2002 Jul 23 (<!> = 4.98). The model UD{A) curves shown as solid, and dotted black 
lines correspond to the identical phases 5.00 and 5.99 from two successive cycles of the the r52(fx) model series, respectively. The model 
spectrum corresponds to the phases 5.00 from the o54(fx) model series. 
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eters are slightly too sraall throughout the 1.1 — 3.8 /im 
region. 

Our data also shows go od agreement with re sults for 
the Mira star S Ori by Wi ttkowski et al.l (|2008[ ). in that 
the minimum UD angular diameters are smallest at 1.3- 
1.4 and 1.6-1.7/J,m and increase by a factor of 1.4-1.5 
around 2.0 /zm. For a full comparison with published in- 
terferomet ric UD angular diamet ers in continuum band- 
passes, see IWoodruff et (l2008h . 

Further work with these data, including asymmetry 
studies, will be presented in subsequent publications. 
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